Low-energy excitations of a linearly Jahn- Teller coupled orbital quintet. 
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The low-energy spectra of the single-mode h®(G+H) linear Jahn- Teller model is studied by means 
of exact diagonalization. Both eigenenergies and photoemission spectral intensities are computed. 
These spectra are useful to understand the vibronic dynamics of icosahedral clusters with partly 
filled orbital quintet molecular shells, for example C60 positive ions. 
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I. INTRODUCTION 

In icosahedral molecules and ions, there occurs the pos- 
sibility of a partial occupancy of a quintet of orbitally 
degenerate electronic levels, of h symmetry. A notable 
example is that of positive Ceo ions, with one (or several) 
holes in a h u molecular orbital 0, • At the adiabatic 
level, partly occupied quintet states are unstable against 
molecular distortions involving the vibrational modes of 
the molecule. Distortions along nondegenerate A modes 
do not involve symmetry reduction, are therefore not of 
Jahn- Teller (JT) type, and can be trivially treated as 
phonon shifts: thus we will ignore A modes altogether. 
Distortions involving fourfold-degenerate G and fivefold- 
degenerate H modes reduce the molecular symmetry, and 
are responsible for the leading (linear) contribution to the 
splitting of the fivefold molecular level. Other (T\ and 
T 2 ) modes are involved in higher-order couplings, which 
become dominating only for substantially large distor- 
tions. Under specific circumstances (null or very strong 
linear coupling) the inclusion of higher-order couplings 
could become essential. However, the present work fol- 
lows the standard practice of neglecting all higher-then- 
linear coupling terms and treating the phonons as har- 
monic oscillators 0, 3 ■ This approximation provides a 
simplified model meant to capture the basic shape of 
the five adiabatic potential energy surfaces (APES) in 
a neighborhood of their high-symmetry conical intersec- 
tion, hopefully extending out to include the JT minima. 
The advantage of the idealized linear model is that it in- 
volves a minimal number of parameters, and that it can 
therefore be studied in detail, drawing fairly general con- 
clusions, applicable to a wide enough range of realistic 
systems. 

Indeed, besides Cgo, orbital quintets could be relevant 
for the spectroscopy of other icosahedral molecules and 
clusters among those synthesized in recent years, includ- 
ing smaller and higher fullcrcncs such as C20 and Cgo, B12 
or B§4 clusters, B12H12, the smallest icosahedral systems 
Sii3, Nai3, or Mgi 3 , and icosahedral structures with 13, 
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19, 55, and 147 atoms in xenon clusters. With different 
coupling patterns possibly realized in those different sys- 
tems, it is of clear interest to investigate how the vibronic 
spectrum changes as a function of electron- vibration cou- 
pling. For simplicity, here we restrict ourselves to a single 
mode (of either G or H type) coupled to the orbital quin- 
tet. 

Several works HH attach the JT problem of the 
orbital quintet in the linear approximation and for a 
single mode. These early works are mainly studies of 
the APES, and in particular of its minima, with the 
phonon coordinates treated as classical variables. Other 
papers deal with the full quantum mechanical prob- 
lem, and provide either approximate |g, |7], |8| or exact 
P, fl(l solutions, sometimes for the many-modes case 
[111 |l2j . The quantum mechanical dynamics of the en- 
tangled electron-phonon system is affected by an elec- 
tronic geometric phase, whose presence/absence has rel- 
evant consequences, in particular, for the ground-state 
(GS) symmetry [1 S E 13 H G3 , thus for magnetic 
resonance spectroscopies [l7| • In particular, Ref. [l(J pro- 
vides a detailed study of the evolution of the tunneling 
splitting through the range of coupling parameters for an 
H mode |3. 

Unfortunately, to date, no experimental determination 
of the tunneling splitting in any icosahedral JT system 
is available. In contrast, spectroscopic techniques can 
and do provide direct measurements of the JT vibronic 
spectrum. For example, molecular photoemission [Tsllill 
was recently shown [ll| to provide relatively detailed in- 
formation on the low-lying part of the vibronic spectrum 
of Cg , especially in a specific symmetry sector [2jj. A 
clean theoretical understanding of the general features of 
this vibronic spectrum is therefore desirable. 

The present work provides precisely an exact determi- 
nation of the low-lying vibronic levels of the single-mode 
linear model, those usually of most direct experimental 
access. The spectrum is studied through a wide range 
of coupling parameters, from weak to strong coupling. 
In practice, the reported results are mainly relevant for 
weak to intermediate coupling, since, as noted above, the 
linear JT model is little more than a mathematical curios- 
ity in the strong-coupling limit, where higher-than- linear 
powers of the distortion become dominant. The vibronic 
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spectra presented in this work illustrate several interest- 
ing quantum phenomena characteristic of dynamical JT. 
Several numerical energies reported in Appendix ^ can 
be used as benchmarks of the accuracy of future calcu- 
lations based on ap prox imate approaches, for example of 
the type of Refs.IlM 

II. MODEL AND CALCULATION 

The model 0, 0, 0] describing the JT coupling of the 
orbitally degenerate h quintet with the molecular vibra- 
tions of symmetry A = G, H is conveniently formulated 
as follows [Hillli: 

H = H o + i?vib + -ffc-v , (1) 
H o = € Y1 ^™<y dm ° ' ( 2 ) 

mcr 

H vih = \t^Y,{ P l + Ql) . ( 3 ) 

i? c _ v = ^ Y, a^ci^n.Q^l^'*- (4) 

ra fimm' 

Here m, label components within the degenerate mul- 
tiplets, for example according to the G5 character in the 
T D D 5 d C5 group chain [Toll23 . G^ 1 , are icosahedral 
Clebsch-Gordan coefficients 24] which couple two ft- ten- 
sor operators and a A tensor operator to a global scalar A 
operator. Q M are the dimensionless normal-mode vibra- 
tional coordinates (in units of the natural length scale of 
the harmonic oscillator), and P M the corresponding con- 
jugate momenta. The multiplicity r = 1,2, needed for 
A = H vibrations only, labels the two separate kinds of 
coupling allowed under the same symmetry |22l l24j : it 
represents the double occurrence of the H representation 
in the direct product h x h (the icosahedral group is not 
simply reducible). Numerical factors k = |5 5 , k H = | 
(which could otherwise be re-absorbed into the definition 
of q\) are included to make contact with the notation of 
Ref.Hi 

Without any loss of generality, the energy position e 
of the h quintet will be taken as the zero of energy. In 
this single-mode problem, all energies are naturally mea- 
sured in units of the vibrational quantum hw. The sec- 
ond characteristic energy scale g\frio (the JT energy) is 
smaller/larger than hu> for small/large g\ (weak/strong 
coupling) . The weak-coupling regime is characterized by 
rapid tunneling among shallow clustered potential wells 
and strong non-adiabaticity associated to the vicinity to 
the conical intersection of the five Born-Oppenheimer po- 
tential sheets; the level structure is perturbatively related 
to the harmonic spectrum of i? v ib- As the coupling grows 
to intermediate and strong, the JT wells deepen and 
move away from the conical intersection and from each 
other: tunneling decreases and the vibronic spectrum ac- 
quires a characteristically intricate structure, which sim- 
plifies again in the strong-coupling limit, where tunneling 



is exponentially suppressed and semiclassical considera- 
tions apply @]. 

The adiabatic approximation, where the vibrational ki- 
netic energy is neglected, provides the main (static JT) 
features of the strong-coupling limit. The distortion op- 
erators Q M are treated as classical coordinates, and op- 
timal distortions are obtained by minimizing the lowest 
APES HH, for example by means of the Opik-Pryce 
method |25j . These optimal distortions (JT wells, or 
minima) are associated to specific splitting patterns of 
the electronic quintet. In particular, G modes produce a 
level pattern of the type (—8, —3, —3, 7, 7), in correspon- 
dence to 10 JT wells of D 3 residual symmetry. 

The situation is slightly more intricate for an H mode 
[HIE EH' due to the multiplicity r = 1, 2. The coupling 
scheme, associated to the r = 2 Clebsch-Gordan coeffi- 
cients, produces the splitting pattern (—4,-1,-1,3,3) 
in correspondence to 10 D3 minima, like the G modes. 
The r = 1 scheme produces a splitting pattern of the type 
(—4, 1, 1, 1, 1), in correspondence to 6 JT minima of D$ 
symmetry. 

In practice, the two-scheme coupling for the H modes 
is a convenient theoretical idealization, since any H mode 
of a real molecule is associated to some amount of cou- 
pling of both r = 1 and r = 2 kind. This admixture is 
conveniently characterized by introducing two coupling 
constants for each H mode, rather than one. The two di- 
mensionless linear coupling parameters gffl can be 
expressed in polar coordinates 0, as 

(!) ( 2 ) • fK\ 

9h -gcosa, gy =g sma. (5) 

In summary, besides its frequency to, each H vibrational 
mode is characterized by a pair (g^\ g^P) of linear cou- 
plings, or equivalently a global coupling intensity g plus 
a mixing angle a. 

Figure shows the evolution of the quintet of elec- 
tronic eigenvalues as a takes its possible values 47] . Here, 
the molecular vibration operators are treated as classical 
variables, and are taken at the relevant static JT mini- 
mum, of -D5 or D% symmetry, appropriate for each value 
of a. Pure coupling of type 1 is recovered at a = and 
7r, while pure coupling of type 2 occurs for a~ir/2. At 
the special angles a s = arctan(3/y / 5) — 0.296 tt ~ 53.3° 
and 7r — a s , where cusps occur in the a-dependence of 
the eigenvalues, D5 and D% minima are exactly degener- 
ate: the JT minima collapse into a single flat trough, and 
the symmetry of the JT system is effectively larger than 
icosahedral 0, 0, IH, HE E3, ■ Figure [H shows that for 
a given coupling energy g 2 huj, the lowering of the lowest 
electronic state (equal twice the static JT energy gain) is 
strongest at a — 0, tt and weakest at a s , ix—a.2- 

To compute the exact vibronic spectrum of H, full 
quantum-mechanical treatment must be applied also to 
the vibrational operators. The standard ladder-operator 
representation for the vibrational coordinates — 

(^bj^ + I \/2 provides a natural product basis of har- 
monic phonon states centered on the undistorted geom- 
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FIG. 1: The splitting pattern of the electronic quintet in 
the h <g> H JT, as a function of the mixing angle a. The 
lowest eigenvalue is nondegenerate, the two upper ones are 
both twofold degenerate; energies are in units of g 2 Kuj. At 
each value of a, Q is taken at the optimal classical distortion 

Qmin • 




JT coupling g G 

FIG. 2: Exact vibronic spectrum of a single G mode coupled 
to a h level, as a function of the linear coupling parameter 
go- (a) GS energy, to which the adiabatic minimum lowering 
Ecins — (2 — j^ga) Fi^> is subtracted, (b) Vibronic excitation 
energies (measured with respect to the GS). All levels are 
labeled according to their symmetry as A (dashed) , T\ + Ti 
(dotted), G (dot-dashed) and H (solid). 



etry, times the electronic quintet. It is straightforward 
to write H matrix elements in this basis and diagonalize 
the resulting matrix. This basis centered on the high- 
symmetry point is especially convenient at weak cou- 
pling, where converged low-energy vibronic eigenenergies 
and eigenstates are obtained by including only few basis 
states. For large g, the number of basis states needed for 
a given accuracy of the low-energy levels grows rather 
quickly, approximately as g 2rl ° a <=, where n osc = 4 for a 
G mode and n osc = 5 for an H mode. For this reason, 
and because of the extreme sparseness of this problem, 
the Lanczos method of diagonalization (2^, H0| is espe- 
cially suitable here. We use a method very similar to 
the one employed in Ref. ITol with the introduction of 
several refinements to deal with the problem of Lanczos 
ghost states, and the symmetry analysis. We check all 
results against basis truncation. For most results, inclu- 
sion of n = 35 phonons in the basis is largely sufficient 
for an accuracy of 10~ 6 huj, but for occasional points in 
parameter space the basis needs to include up to n = 60 
phonons (over 8 million states) to guarantee the required 
accuracy. The spectra computed with this method and 
reported below draw a full quantitative link from weak 
through intermediate to strong coupling for the h ® G 
and h® H linear dynamical JT model. 



III. VIBRONIC SPECTRA: ENERGIES 
A. The h ® G JT 

Figure [2^. shows the difference between the exact dy- 
namical JT GS energy and £ c ias (the classical adiabatic 
energy —j^g^ku) of the 10 D3 wells plus the zero-point 
energy Ihio of the uncoupled oscillators), as a function 
of the coupling strength go to a G mode. This residual 
energy represents the part of JT energy gain not ascrib- 
able to trivial adiabatic lowering of the JT minima with 
respect to zero coupling: this energy difference gauges 
quantum-mechanical effects on the vibrational motion. 
The rapid decrease in GS energy at weak coupling is due 
to the large drop in quantum kinetic energy in the transi- 
tion from harmonic oscillations around one isolated adi- 
abatic minimum to rapid tunneling among 10 shallow 
minima. As these minima move apart, tunneling is sup- 
pressed, and the plotted energy difference turns back up 
to the strong-coupling value (3" 1 / 2 — l)Hu>, which mea- 
sures the zero-point energy gain due to phonon softening 
at the anisotropic JT minima Indeed, the second- 

order expansion of the lowest APES around any of the 
equivalent JT minima yields two normal modes (of sym- 
metry Ai and Ai in the local D3 group) both of frequency 
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u>, plus a twofold-degenerate softer mode (of symmetry 
E) of frequency 3 -1 / 2 tJ. The difference in harmonic zero- 
point energy 2 \huj + 2\2>- 1 / 2 huo-A\hui = {i" 1 / 2 - l)hw 
provides an accurate estimate for the exact vibronic en- 
ergy (Fig. at strong coupling, where inter-well tun- 
neling is suppressed and the phonon dynamics reduces 
to harmonic oscillations around the JT wells. Note that 
the "new" frequencies at the wells are independent of the 
coupling strength go, since for different coupling go, the 
set of five APES differs only for a scaling factor, thus 
remaining self-similar. At weak coupling, the harmonic 
frequencies at the JT wells remain unchanged, but this 
expansion is practically irrelevant, since quantum kinetic 
energy delocalizes the motion far enough from the min- 
ima for anharmonic and nonadiabatic effects to domi- 
nate. 

The vibronic GS symmetry remains H at all couplings. 
This is eventually a consequence of a Berry phase of n 
acquired by the electronic wavefunction as the distor- 
tions follow any of the cheapest (5-corner) loops through 
adjacent -D3 minima: this entanglement makes the to- 
tally symmetric A vibronic combination of the wells less 
advantageous than the H combination, which thus re- 
mains the GS for all go- Indeed, as seen in Fig. Eb, even 
at strong coupling, the A vibronic state lies above both 
the H GS and a G excitation At strong coupling, these 
5(H) + 4(G) + 1(A) states represent the 10 symmetry- 
adapted vibronic combinations of the harmonic GS's at 
the 10 D3 wells. Likewise, in the strong-coupling limit, 
higher excited states, represent suitably symmetrized 
combinations of vibrational excitations in the wells. In 
particular, the 20 states (Ti + T 2 + G + 2 H) converging 
at strong coupling to 3 _1 / 2 ?iw represent the tunneling- 
symmetrized one-phonon oscillations of the new softer 
frequency. Similarly, the 20 states (A + T\ + T 2 + 
2 G + H) converging at strong coupling to fkv represent 
the tunneling-symmetrized one-phonon oscillations of the 
modes at the original frequency hio. 

At weak coupling, the vibronic levels reconstruct the 
harmonic ladder, as expected when H c _ v is a weak per- 
turbation. The present exact calculation traces precisely 
the crossover from the weak to strong coupling regime: 
for accurate numeric eigenvalues, see Table [i] in the Ap- 
pendix. It is worth noting that the lowest A level, which 
at strong coupling is a low-energy tunneling excitation, 
correlates to the 2 ha; multiplet at weak coupling Q • Fi- 
nally, note that the Ti and T 2 levels are exactly degener- 
ate for all values of the coupling. 



B. The h ® H JT 

We move on now to study the spectrum for a vibra- 
tional mode of H symmetry. The Hamiltonian H de- 
pends now on two parameters. Each vibronic eigenvalue 
traces a 2-dimensional surface as g and a (or equivalcntly 

Oh and gjp) are varied. The most effective presentation 
of these spectra is obtained by means of radial cuts at 
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FIG. 3: "Tangential" normal frequencies ui T of vibration 
around the JT wells, in units of the original frequency ui, 
as a function of a. Both these modes are twofold degenerate. 
As customary of linear JT, the frequency of the fifth (radial) 
normal mode equals to, regardless of a. 

fixed a, plus circular cuts at fixed g. 

Like for the G mode, the strong-coupling asymptotic 
vibronic energies, are determined by the normal frequen- 
cies of oscillation around the adiabatic JT minima, and 
these frequencies now depend on a. The five normal 
modes separate into a radial vibration of unchanged fre- 
quency u> and symmetry Ai, plus two softer doubly- 
degenerate (of E\ and E2 symmetry for D5 wells and 
both of E symmetry for D3 wells) tangential modes, rep- 
resented in Fig. |3| For special a values, extra degenera- 
cies occur: 

• for a = and 7r, the two tangential modes become 
degenerate at lo t — (9 / 40) 1 1 ' 2 oj; 

• for a ~ 0.1347T and 0.866 7r, the upper tangential 
frequency tops at a maximum where it reaches the 
frequency uj of the radial mode; 

• at the special a = a s and tt — a s points, both 
tangential frequencies vanish as the vibrations turn 
into free modes of pseudorotation along the flat JT 
trough; 

• in addition, for a = 7r/2, both soft modes reach a 
local maximum, of frequencies u>r — (5 /lA) 1 / 2 ^ 
and u;i 2) = (5/6) ^w. 

Observe that the normal frequencies at the wells only 
depend on a, while they are independent of the coupling 
strength g, like for the G mode. Again, these single- 
well harmonic frequencies determine the strong-coupling 
asymptotic vibronic energies, but they are essentially 



5 




JT coupling g 



FIG. 4: Exact vibronic spectrum of a single H mode coupled 
to a h level, as a function of the linear coupling parameter 
g, for a — 0. (a) GS energy, to which the energy of the 
adiabatic JT wells i? c i as = (5/2 — g 2 /10) Hlu is subtracted, 
(b) Vibronic excitation energies (measured with respect to 
the GS). The levels are labeled according to their icosahedral 
symmetry even though, for this special value of a, the effective 
underlying symmetry is spherical. The conversion in spherical 
angular- momentum labels [j] is A = [0], H = [2], G + H = [3], 
T1+T2+G = [1] + [3]. (c) Comparison of the exact excitation 
energy of the lowest A state (the tunnelingsplitting) with the 
approximate expression obtained in Ref. by means of the 
unitary transformation method [32^|. 



irrelevant at weak coupling, where inter-well tunneling 
dominates. 

Not unlike the h®G calculation, to draw the GS energy 
it is convenient to subtract the quantity 

![| — jij(coscv g) 2 ] fad, < a < a s 
■ ( 6 ) 
[§ _ i^(sina g) 2 ] fau, a s <a<\ 

E c i as (g,a) represents the trivial zero-point energy \fao 
of -f/vib plus the adiabatic lowering of the D 5 / D 3 JT well 
appropriate to that value of a. 

Figure 0Ji reports the exact GS energy, for the spe- 
cial mixing angle a = (a mode of pure type r = 1). 
This energy difference converges rapidly to its asymp- 
totic value [(18/5) 1 / 2 — l}fao at strong-coupling. Like for 
the G mode, this is understood in terms of zero-point en- 
ergy associated to the four frequencies uj t — (9/10) 1 ^ 2 uj 
plus one frequency to of the normal-mode vibrations at 
the JT minima. 



The spectrum of excitations Fig. 0Jd is much less clut- 
tered than that of the G mode (Fig. I^Jd). The main 
reason is the larger than icosahedral effective symmetry 
of the Hamiltonian at this specific coupling angle. This 
larger (spherical) symmetry induces extra degeneracies 
among the vibronic levels, which one could in fact label 
with spherical angular- momentum labels. The lowest ex- 
cited state is of A symmetry: it is the one state which 
at strong coupling drops towards the H GS. Together, 
these 5(H) + 1(A) states represent the 6 symmetrized 
tunncling-split vibronic combinations of the harmonic 
GS's at the 6 -D5 wells. An approximate analytical ap- 
proach to this problem based on this idea (the unitary 
transformation method [22) yields an estimate of 
the excitation energy (solid line in Fig.^Ji) 4{|: this ap- 
proximation is seen to be fairly accurate throughout the 
range of couplings. The approximate expression underes- 
timates the tunneling splitting at strong coupling due to 
the simplifying assumption of isotropic JT wells, of fre- 
quencies all equaling lj, which therefore exaggerate the 
localization of the distorted states. Corrections due to 
the anisotropy of the wells were introduced in Ref. I2H 
In the Appendix, Table [H] reports a few of the vibronic 
eigenenergies used to draw Fig. 0] 

Higher excitations cluster around the harmonic 
normal-mode energies at the D5 wells: the tangential 
fao T = (9/10) 1 ^ 2 fav and the radial fao. The 24 states 
(Xi +T2 + 2G + 2H) around fao T are symmetrized com- 
binations of the four 1-phonon states per each of the six 
D5 JT wells. Likewise, the 6 states (H + A) near fao 
are combinations of the 1-radial-phonon states in the 6 
wells. Note in particular that the two 1-phonon H states 
remain almost degenerate until g < 2, they split signif- 
icantly in the interval 2 < g < 8, and then re-converge 
again (to fao T ) at strong coupling. At strong coupling, 
higher excitations (not drawn) are found in the overtone 
region 2 fav T , and above. 

For a = 7r/10, Fig.OJ, shows the exact GS energy re- 
ferred to E c \ as . This small value of a has little effect on 
the GS energy, as compared to the a = case. On the 
contrary, the comparison of Fig. \^jp with Fig. shows 
that even such a small value of a (i.e. a tiny admixture 
of the coupling of type r = 2) is sufficient to completely 
resolve the "accidental" degeneracies of the a — spec- 
trum. A and G states are little affected, while the Ti/T 2 
degeneracy is widely broken. Likewise, the weak-coupling 
near degeneracy of the 1-phonon H states for a = is 
now widely split. Eventually, at strong coupling, these 
two H states converge to different tangential frequencies 
of the JT wells. 

Coming to the spectrum for a = 7r/5, Fig. shows 
a continuous evolution away from a = through a — 
7r/10. We observe that, at strong coupling for both a = 
7r/10 and a = 7r/5, only one new frequency is seen in 
the strong-coupling spectrum: that of decreasing uj t for 
increasing a (u) T — 0.86 uj for a = tv/10, lu t = 0.69 a; 
for a — 7r/5). Figure |21 confirms that the frequency uj t of 
the other tangential mode is very close to its maximum u>, 
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FIG. 5: Exact vibronic spectrum of a single H mode coupled 
to a ft level, as a function of the coupling parameter g, for 
a — 7r/10. (a) GS energy, to which -E c i as is subtracted, (b) 
Vibronic excitation energies measured with respect to the H 
GS. All levels are labeled according to their symmetry as A 
(dashed), 7\ / 2 (dotted), G (dot-dashed) and H (solid). 
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FIG. 6: Same as Fig. for a — tt/5. 
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FIG. 7: Same as Fig. 0for a = 0.28 tt < a 3 



thus very difficult to separate on this scale. In practice all 
18 relevant one-phonon states converge near Hu> at strong 
coupling. 

Moving further away from this accidental degeneracy, 
to a — 0.28 7T (Fig. 0|, the two new tangential frequen- 
cies become uj^ — 0.33 lu and oj^ ~ 0.61 u. These 
frequencies are so small because of the nearness to the 
special angle a s where -D5 and D3 minima degenerate to 
a flat trough (see Fig.[3J). These small values explain the 
twofold reason why the spectrum of Fig. 0) is so clut- 
tered: (i) tunneling is only weakly suppressed by the low 
inter-well barriers, and (ii) several overtones and combi- 
nations of these soft tangential modes fit in the consid- 
ered energy range. The large zero-point energy lowering 
shown in Fig. 0l is also well accounted for by the values 
of u> T . 

The next point a = 37r/10 (Fig. [SJ is located imme- 
diately beyond a s , with D3 minima barely prevailing 
over D§ distortions. Here, the tangential frequencies are 
as low as ~ 0.125 uj and uj^ — 0.247 u>: this ac- 
counts for the huge zero-point energy lowering apparent 
in Fig. [SJl in the strong-coupling limit. In contrast to 
all previous a values, where D5 wells prevail, we observe 
here the celebrated jf| 0, EH level crossing (occurring 
at a rather strong coupling g ~ 12) to a nondegenerate 
GS of A symmetry. The vibronic spectrum (Fig. [Bp) is 
now even more cluttered than for a — 0.28 n, it shows 
an intricate pattern of level crossings (states of different 
symmetries) and avoided crossings (states of the same 
symmetry). Several higher- lying levels have been omit- 
ted in the plot for clarity. This is due to the exceedingly 
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FIG. 8: Exact vibronic spectrum of a single H mode coupled 
to a h level, as a function of g, for a = 0.3 n > a s . (a) GS 
energy, to which -E c i as is subtracted; the GS is of H symmetry 
at weak coupling, and of A symmetry at strong coupling, 
(b) Vibronic excitation energies measured with respect to the 
lowest H state. 



low normal frequencies at the wells, plus the larger num- 
ber of wells (10 instead of 6). The clustering around the 
semiclassical frequencies is barely hinted, even for the 
rather strong coupling at the right side of the plot. 

Coming now to a = 2ir/5 (Fig.^), the tangential fre- 
quencies at the -D3 minima are now u>t — 0.529 u and 

(2) 

lut ~ 0.854 u>, much larger than for the previous cut 
at a = 37r/10. These values account for the reduced 
zero-point energy lowering shown in Fig. Also, the 
spectrum (Fig. ^jp) is less dense: the 10 tunneling states 
(A + H + G) converge to zero energy at strong coupling; 

the 20 1-phonon tangential states converging to fouir are 

(2) 

fairly well visible; the tangential hw r , radial hw, and 
overtone 2 ftw^ states are still substantially intermixed 
even at the large coupling of the right side of Fig. |5Jd. 
The splittings between T\ and T 2 pairs are consistently 
reduced here, compared to previous values of a. Note in 
particular that, in contrast to the a range producing D§ 
wells, both 1-phonon H states correlate at strong cou- 
pling to the lowest tangential mode of frequency hu>^\ 
A fairly general feature of the vibronic spectrum, espe- 
cially visible in Figs. |2t>, Eb ; and ITUb . is the rapid drop 
of several vibronic energies, then followed by a successive 
climb back at strong coupling where tunneling is sup- 
pressed and the motion gradually collapses to oscillations 
around the adiabatic wells. 



FIG. 9: Same as Fig. 03 with a = 2tt/5. 

At a = 7r/2 (Fig. EH, JT coupling of pure r = 2 type 
is acting. As illustrated in Fig. |3 above, the tangen- 
tial frequencies at the JT wells reach here a local max- 
imum wf' ~ 0.598 uj and ~ 0.913 cj. Accordingly, 
the zero-point energy lowering of Fig. HOfa , amounts to 
[(5/14) 1 / 2 + (5/6) 1 / 2 - 2}hw = -0.490 hw. The GS level 
crossing occurs now at a marginally smaller g ~ 10.7, 
due to suppressed tunneling splitting. The spectrum is 
now even less dense than for previous values of a, due 
to the larger frequencies at the JT wells, and the exact 
degeneracy of the T\ and Ti pairs. The 20 1-phonon tan- 
gential states (G + 2H + T\ + T%) converging to hu>^ are 
well visible; the 20 states (of the same symmetries) con- 

(2) 

verging to hiOr are still entangled with the 10 huj radial 
states (A + H + T 1 +T 2 ). 

For 7r/2 < a < tt the spectra are mirror symmetric 
around a = ir/2 to those reported for the interval < 
a < ir/2. The only difference is that all T\ and Ti states 
(which always come in pairs) are exchanged. 

The cv-dependence of the low-lying excitations (mea- 
sured with respect to -Egs) is best illustrated by fixed-g 
circular cuts. Figure UTI reports the exact excitation en- 
ergies at fixed g = 1 for the 25 vibronic states related to 
the 1-phonon excitation and for the 75 states related to 
the 2-phonon overtone, as a function of a. This value of 
g is sufficiently small for perturbation theory to provide 
not especially bad estimates of the excitation energies 
drawn in Fig. 1111 A, G, and H energies are symmetric 
around a = tt/2, while T\ and T 2 energies are antisym- 
metric functions of a, coming in pairs. The lowest A 
level is the lowest excitation for all a: at strong coupling 
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10 

JT coupling g 

FIG. 10: Same as Fig. |U with a = tt/2. 



15 




FIG. 11: Exact vibronic spectrum of a single H mode cou- 
pled to a ft electronic quintet, as a function of a, for g = 1. 
(a) The vibronic states related to the first overtone 2hui re- 
gion, (b) The vibronic states perturbatively related to the 
"fundamental" excitation at hui. 




FIG. 12: Low-energy region of the vibronic spectrum of a 
single H mode coupled to a ft level, as a function of a, for 
g — 11. Energies are referred to the lowest H state. 



it evolves into a low-energy tunneling partner with the 
H state (and a G state when D% minima prevail). The 
average position of the two H states in the one-phonon 
multiplet is located above Huj: this causes the characteris- 
tic signature of JT in molecular photoemission which was 
observed and discussed in Ref. [llj. Note that no singular 
behavior is observed near a s and it— a s (vertical lines), in 
contrast to all quantities based on a classical treatment 
of the vibrations, such as those plotted in Figs. ^ and [3] 
the transition from D5 to D3 wells is perfectly smooth in 
the quantum spectrum, and signaled by only a number 
of extra degeneracies. 

Figure^]reports the excitation energies of several low- 
lying vibronic states for g = 11 (intermediate/strong cou- 
pling). The tunneling partners are visible at the low- 
energy end, with the crossing to an A GS in the range 
0.427T < a < 0.587T. The coupling is not strong enough to 
have the higher vibronic levels organized according to the 
harmonic frequencies of Fig. [21 throughout the a range, 
but similar trends are already recognizable: in particu- 
lar for a < 0.27T and a > 0.87T, the vibronic levels follow 
quite closely the harmonic frequencies of Fig.Ql and their 
overtones and combinations. As observed in the discus- 
sion of Fig. El the vibronic spectrum in the strong tunnel- 
ing regions at and around the flat-trough values a = a s 
and 7T — a s is extremely congested: only the lowest-lying 
states are drawn in figure. Accordingly, most vibronic 
levels undergo a fast movement in that region, but al- 



9 



ways evolve continuously across a — a s and n — a s . In 
the limit of infinitely large g, the singular behavior of 
Fig. would eventually emerge even in the exact quan- 
tum spectrum. The very large slope in the a dependence 
of many vibronic energies may prove useful in future pre- 
cise determinations of the JT coupling parameters from 
spectroscopic data. 



IV. PHOTOEMISSION SPECTRA 

Several spectral properties of experimental interest are 
readily accessed by the exact diagonalization method at 
hand. Photoemission spectra (PES) involving the quintet 
h molecular orbital are of course affected by JT coupling 
to the degenerate vibrations. If, as is usually the case, the 
photoemitted electron kinetic energy is much larger than 
ftu>, photoemission occurs in a time much shorter than 
that characteristic of phonon dynamics. In this limit, it 
is convenient to apply the sudden approximation |.'i.'ll3 1| . 
where the photoemission process is described simply by 
the operator c ma suddenly destroying a spin-u electron 
in orbital component m of the quintet level. The phonon 
shakeup contributions to the PES are obtained in terms 
of the matrix elements of this hole-creation operator c m „ 
between the initial configuration and the final vibronic 
states. If we assume that the initial temperature is neg- 
ligible (fceT <C h\o, only the GS |GS) is initially popu- 
lated) , then the PES intensity for the creation of a spin-cr 
hole in orbital m is 

I ma (E) cx^|(/|c mCT |GS)| 2 6(E-E f ) , (7) 
/ 

where averaging over the components is implied when- 
ever |GS) is degenerate. If the individual contributions 
of different orbital and spin component are not separate, 
the total PES intensity is then 
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FIG. 13: Zero-temperature PES obtained by fast removal 
of the single electron from a h orbital coupled to a single H 
mode, for several values of the coupling g, and a = it/5. En- 
ergies are referred to the GS; a phenomenological broadening 
HWHM=0.02 hui provides a finite width to the peaks. Frac- 
tional spectral weights of the 0-phonon and 1-phonon line are 
listed. 



all non-adiabatic effects and phenomena neglected in the 
Frank-Condon approximation |36| are fully included. 



JT in the initial states 



For generic occupancy of the h level, the spectrum 
is generally affected not only by JT, but also by intra- 
shell electron-electron repulsion, and the ensuing multi- 
plet structure 0, 12^1 l35l | : calculation of the PES in such 
conditions must take into account several other system- 
specific interactions. It thus goes beyond the study of the 
single-mode rather general model considered in this work. 
However, electron-electron repulsion plays no role in two 
relevant cases, which are therefore completely described 
by the JT Hamiltonian H of Eqs. Q-lQ}: (i) when a 
single electron occupies the quintet orbital in the initial 
state, and (ii) when the orbital quintet is initially com- 
pletely full and photoemission produces a single hole in 
the final state. In the first case, JT only affects the initial 
configuration, while in the second case, JT only affects 
the final states. The application of exact diagonaliza- 
tion to the two cases separately provides spectra where 



When JT only affects the initial states, the molecule is 
instantaneously brought from the multi-well intrinsically 
nonadiabatic vibronic GS to a linear superposition of the 
final harmonic oscillator states. According to Eqs. Q 
and JSJ), the peaks in the spectrum are located at the 
final-state energies Ef, which here simply involve sin- 
gle or multiple excitations of the harmonic phonon fuj. 
The spectrum is therefore a sequence of regularly spaced 
peaks, with intensity of the 0, 1, 2... -phonon line propor- 
tional to the square modulus of the corresponding total 
component in the initial vibronic GS. 

Figure presents precisely a spectrum with this sim- 
ple structure, for an H mode characterized by a = 
7r/5. According to the standard electron-phonon (non- 
JT) Frank-Condon picture, the percentage of spectral 
intensity displaced away from the 0-phonon to the multi- 
phonon excitations increases with the distance between 
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FIG. 14: Same as Fig. HSlfor a = 2vr/5. 



the initial and final equilibrium configurations, which, in 
turn, is proportional to the linear coupling parameter g. 

The situation is qualitatively similar for a = 2n/5 (and 
for the case of a G mode, which is omitted for brevity). 
The only quantitative difference to observe in Fig. ITU 
compared to Fig. Ellis that for a = 2tt/5 spectral weight 
is preferentially transfered to the 1-phonon line, while 
for a — 7r/5 it is more equally distributed among several 
phonon excitations. 

The spectra in Figs. 1131141 are obtained by means 
of a complete Lanczos calculation of the GS wavefunc- 
tion. This is straightforward, as long as zero temper- 
ature is addressed as in the present work. If instead 
finite-temperature was to be considered, many initial 
vibronic states would be required for a thermal aver- 
age, and this is generally a difficult task for the Lanczos 
method [37]]: completely different techniques, e.g. based 
on Monte Carlo, could prove more convenient there. 



B. JT in the final states 

When JT affects the final states, the energies Ef mark- 
ing the PES peaks are the vibronic energies of Figs.|21E} 
1121 The transparent structure of a harmonic spectrum is 
now lost. However, the calculation of the PES is straight- 
forward, as |GS) is now the 0-phonon state of the fillcd- 
shell non-JT molecule The symmetry of the state 
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FIG. 15: Zero-temperature PES obtained by fast removal of 
an electron from a filled h orbital coupled to a single G mode, 
for several values of the coupling g. Energies are referred to 
the GS; a phenomenological broadening HWHM=0.01 huj pro- 
vides a finite width to the peaks. Fractional spectral weights 
of the "zero-phonon" line are listed. 



c-mcr | GS) occurring in the matrix element of Eq. Q is that 
of the electronic operator c mcr , namely H. As a conse- 
quence, only transitions to H final states (/| are possible, 
while matrix elements to all other vibronic states vanish 
by symmetry. Peaks occur therefore only at the energies 
marked by the solid lines in Figs. 121 141121 

Figure ITS1 reports the computed PES for a coupled G 
mode and a few intermediate values of g. The positions 
of the low-lying peaks correspond to the solid lines of 
Fig. [3d, and are readily followed. Strong intensity trans- 
fers due to level crossings are well visible. The intensity of 
the "zero-phonon" line (better referred to as the vibronic 
GS) decreases steadily as g increases [50|: the intensity 
lost there redistributes across an intricate vibronic spec- 
trum. 

Figures ITTdI and IT7I report the computed PES for a cou- 
pled H mode characterized by the same value of a of 
Fig- ED producing D$ minima) and of Fig. |SJ (2tt/5, 
producing D3 minima) respectively. The positions and 
intensities of several low-lying peaks can be traced: sev- 
eral intensity transfers due to level crossings are well vis- 
ible. Again, the intensity of the "zero-phonon" line de- 
creases consistently as g is increased, while the intensity 
lost there redistributes across a rather dense spectrum of 
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FIG. 16: Zero-temperature PES obtained by fast removal 
of an electron from a filled h orbital coupled to a single H 
mode, for several values of the coupling g and, a = n/5. En- 
ergies are referred to the GS; a phenomenological broadening 
HWHM=0.01 huj provides a finite width to the peaks. Frac- 
tional spectral weights of the "zero-phonon" line are listed. 



vibronic levels. 

Figures 1151171 illustrate in detail the general finding 
[TT| that in the PES of dynamic JT system with weak to 
intermediate coupling, the "1-phonon line" splits into es- 
sentially two features, usually with most spectral weight 
in the upper peak, blue-shifted up to 30% higher energy 
than Huj. The blue shift of the 1-phonon satellite is not 
specific to the h<g> (G + H) JT model. For example, sim- 
ilar conclusions apply to the e®£ problem [33 ■ Indeed 
a 15% blue shift of the first vibronic satellite is clearly 
observable in the PE spectra of benzene (3^ • 

A general feature of the displayed PES is that the per- 
turbative structure of n-phonon states, still fairly well 
recognizable for g — 2), is rapidly lost in the PES for 
g > 3. The only other general consideration we can draw 
from these spectra is that not much regularity is to be 
expected: unless coupling is very weak, phonon shake- 
ups in the PES from a JT molecule build very intricate 
spectral structures, which are strongly g-dependent (also 
a-dependent, for H modes), even for a single phonon 
mode. 
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FIG. 17: Same as Fig.EOfor a = 2tt/5. 



V. DISCUSSION 

Understanding the general features of a vibronic spec- 
trum is a necessary step towards the interpretation and 
assignment of the spectra of molecules and molecular 
ions. This has been undertaken and carried out quite 
extensively for electronic doublets, triplets ]38l lifl Jill 
and more complex electronic configurations |35l |42j, |43j ■ 
The present paper fills a gap in the theory of icosahe- 
dral systems, providing general guidelines for the inter- 
pretation of vibronic PES of the electronic quintet in the 
h ® (G + H), restricting to the idealization of a single 
mode. In practical applications, the situation is much 
more intricate. All actual icosahedral systems are many- 
mode systems, where the inter-mode interaction gener- 
ates a much less coherent spectrum than simple non-JT 
phonons |44|. The present theory allows the basic un- 
derstanding of the PES of icosahedral systems wherever 
the coupling to a single high-energy mode is dominant, 
the other couplings only contributing to the background. 
The couplings of fullerene Cgo are approximately of this 
kind, and the spectral of Sect. IIVBI allow to rational- 
ize some observed PES features in terms of an effective 
single-mode theory. 

On the other hand, the spectra of Sect. II V XI describe a 
situation where a single electron in a h state is rapidly re- 
moved. This cannot occur in practically accessible ionic 
states of Ceo, but could well occur in singly charged neg- 
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ative ions of higher icosahedral fullerenes, under the con- 
dition that the LUMO is a quintet. By particle- hole sym- 
metry, the same theory applies to inverse PES, where a 
single hole in a quintet is filled. This is relevant, e.g. for 
inverse PES of C£ . 

The difficulty of the many-modes JT [H HI] is that 
no linear superposition of the spectra of single modes oc- 
curs (except for the weak-coupling limit). This means 
that any quantitative description of the dynamics re- 
quires to include all modes together in a single calcula- 
tion. The method of Lanczos exact diagonalization used 
in the present work is b y n o means limited to the single- 
mode case. Indeed, Ref . Hll deals successfully with the in- 
tricacies of the Ceo many-mode problem, also accounting 
for thermal effects. The same theory used in Sect. IIV Bl 
for PES could be applied to an optical transition from 
a nondegenerate electronic a state to a final h state, if 
that dipole-forbidden transition could be realized. Modi- 
fications of the same method could be employed to study 
phonon shakeups in optical absorption/emission in dipolc 
allowed transitions involving electronic degenerate states 
at both sides: these are often associated to "product" 
JT systems [3^, , and with due account of "term" ex- 
change interactions, JT would account for the vibronic 
"decorations" of the spectra, not unlike PES. More ap- 
plications of the Lanczos method (including generaliza- 



tions thereof |3(J) to other spectroscopical applications 
are therefore expected in the future. 
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APPENDIX A: EXACT EIGENENERGIES 

Tables Ull VI collect a number of exact vibronic eigenen- 
ergies of the orbital-quintet single- mode linear JT model. 
Note that .Ecias was not subtracted here, therefore the 
tabulated energies represent precisely the eigenvalues of 
the Hamiltonian as defined in Eq. (JIJ. All reported dig- 
its are significant. Extra degeneracies are apparent in all 
tables except Table ITTT1 (a — 7r/4). However, the appar- 
ent degeneracy between the A states and H states for 
g = 15, a = (Table [HJ only indicates that the (tun- 
neling) splitting between these states is smaller than the 
energy resolution of the diagonalization (see also Fig.0J. 
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TABLE II: Exact eigenenergies of several low-lying vibronic 

states of h <g) H linear JT for a — and selected values of the TABLE IV' Same as Table ITTl for a = tt/2 

coupling intensity g. Energies are in units of hu). 
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